Inclusive hadron production at LHC in dipole momentum space 
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Within the fct-factorization formahsm and using a momentum space model for the dipole scattering 
amplitude we describe the LHC data on proton-proton collisions, and also on the ALICE preliminary 
data on the 5.02 TeV proton-lead run. The model is based on the asymptotic solutions of the 
Balitsky-Kovchegov equation, being ideal in the saturation domain where the target wavefunction 
has high occupation number. We compare the results of this work with the analysis performed 
recently in the scope of a hybrid factorization formalism. 



I. INTRODUCTION 

The ALIGE collaboration at LHC has just released [ij 
preliminary data on proton-nucleus collisions, which will 
be crucial on the understanding of the hadronic matter 
forming the Quark Gluon Plasma in the case of nucleus- 
nucleus collisions. This is so because on the p -(- Ph run 
the presence of an intricate final state is not expected, 
being this sort of collision important to study the initial 
state phenomena and separate these from those intrinsic 
to the final state. In such a high energy collision, the 
initial state is expected to be a many body system whose 
dynamics can be described through the Color Glass Con- 
densate (CGC) formalism (for reviews see 0-13 )• In this 
framework, gluon recombination processes should hap- 
pen in addition to radiation one, in order to saturate the 
fast growing of the cross section and thus keep the scat- 
tering amplitude unitary. The main aspect of such sat- 
uration formalism is the presence of a dynamically gen- 
erated scale, the saturation scale Qs{Y), marking when 
the recombination effects start to become relevant during 
the evolution on the rapidity variable Y = ln(l/a;) [11-01 ■ 

The simplest way to look for such phenomena is us- 
ing the Balitsky-Kovchegov (BK) evolution equation for 
the color dipole scattering amplitude 0. It resumms 
all the powers of Q!slog(l/x) in the LLx approxima- 
tion, and includes a nonlinear term responsible for the 
gluon merging in the target wave function. By the mo- 
ment the numerical solutions of BK equation are known 
at NLO accuracy, giving rise to a numerical model of 
the dipole scattering amplitude Q. Here, however, we 
will focus on the phenomenology of its asymptotic solu- 
tions, which are easily obtained through a mapping of 
the BK equation into the Fisher-Kolmogorov-Piscounov- 
Petrovsky (FKPP) [l^l reaction-diffusion equation. This 



statistical formulation allowed us to see the universal be- 
havior of the high energy QCD amplitudes. One knows 
that the saturation scale increases with a power law of 
the rapidity in the LO, while behaves as ^/Y in the NLO 
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Other striking property of such formulation 
is related to t he g eometrical scaling property of the DIS 
cross sections [ijl . Such property, observed in the HERA 
DIS data, establishes that in the very high energy limit 
the cross section depends on only one variable, the ra- 
tio Q^/Q^(F), instead of both and Y variables sep- 
arately, making clear the importance of the saturation 
scale on the observables at high energies. It turns out 
that such scaling behavior can be seen as the formation 
of a traveling wave pattern in the solutions of the BK 
equation, in the analogy with the FKPP equation [lTj |. 

Some phenomenology has been done in the last years 
using the traveling wave description of the high energy 
QCD amplitudes, but here we will focus on the AGBS 
model for the dipole amplitude [l^. This model was 
already used to study possible presence of the fluctua- 
tion effects on the gluon number in DIS data [l3|, as 
well as to investigate the nuclear effects present on the 
small- a; region [1^. Recently it has been shown that 
the AGBS model describes equally well the DIS and in- 
clusive hadron production in p -f ^ and p -\- p collisions 
[l6| . This was made in a new simultaneous fit to HERA 
data on the proton structure function [T^ and RHIC data 
on the pt distribution of the produced hadrons and pi- 
ons from BRAHMS and STAR collaborations [H 
The resulting fit, which is based on the hybrid formal- 
ism and uses the parton distribution functions (PDFs) 
for the hadronic content of the projectile, also describes 
equally well the LHC data on the hadron yield for p+p 
collisions, although with large K factors. Such formal- 
ism is good to describe the proton fragmentation region 
in the forward direction, where the nucleus is probed at 
very small-x while the proton's values of x are larger, 
however it is not the best choice in the central rapidity 
region. Here both hadrons wave functions are probed 
at small- a; , so that one can employ the fct-factorization 
formalism where the colliding hadrons are described in 
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terms of their unintegrated gluon distribution (UGDs), 
or simply unintegrated dipole distribution in the case of 
the BK evolution for the scattering amplitudes. 

In this work we confront the AGBS dipole model with 
the LHC data on single inclusive charged hadron produc- 
tion from p + p reactions using the fct-factorization for- 
malism. We also apply the model to deutcron(proton)- 
nucleus collisions at RHIC (LHC), comparing the predic- 
tions with the first data on p + Pb measured recently by 
the ALICE collaboration. 

The paper is organized as follow: in section|II]we revise 
the LO BK equation in momentum space and its traveling 
waves solutions giving rise to the AGBS model. Section 
mil is devoted to establish the formulae and physics be- 
hind the fct-factorization formalism we use in this work. 
In section HV] we show the results and a discussion is left 
for section Ivl 



II. BK SOLUTIONS AND THE AGBS MODEL 

The Balitsky-Kovchegov (BK) equation Q is a non- 
linear evolution equation in the rapidity variable Y = 
ln(l/a;) for the forward amplitude N{r, Y) of a qq dipole 
scattering with coordinates r = x — y oE a. target. It 
can be derived on the Mueller's dipole picture |2(| for 
the high energy scattering, that uses the large limit 
approximation proposed by t'Hooft in the seventies and 
thus guaranty that the evolved dipoles interact indepen- 
dently with the target. When neglecting the impact pa- 
rameter dependence, one can write the dipole amplitude 
in momentum space as 



(1) 



with r being the transverse size of the interacting dipole, 
while k denotes its reciprocal transverse momentum. The 
BK evolution equation for this amplitude reads 



dYN{Y, k) = ax{-dL)N{Y, k) - aN{Y, kf, (2) 



where 



X(7) = 2^(1) - ^-(7) - V'll - 7) 



(3) 



is the BFKL [26| kernel defined in terms of the digamma 
functions, iplj) = r'(7)/r(7), while L = log(fcVfcg), be- 
ing ko a fixed soft scale. 

The asymptotic solutions of such equation can be 
obtained through a map of the BK equation into 
the Fisher-Kolmogorov-Piscounov-Petrovsky reaction- 
diffusion equation [lO| , known to admit traveling waves as 
solutions. This imply that the wavefronts have the form 
f{x — Vet) for large values of fc, which translates into the 



geometrical scaling property of the BK amplitudes 



N {Y, k) 



log 



Xoj {k^/Ql{Y)) 
2ax"{lc)Y 



X exp 



(4) 



where the saturation scale is defined as 



\{Y)^Qlc^^[\Y 
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with A = OLx{lc)Hc measuring how fast the amplitudes 
reach the saturated domain. 

Within this traveling wave description of the high en- 
ergy QCD amplitudes, a model for the dipole scatter- 
ing amplitude was proposed, which describes the in- 
frared behavior as a Fourier transform of a step function 
QirQs — 1), demanding unitary amplitudes in the coordi- 
nate space. Thus, in this region the amplitude is written 
as 



N{Y,k) w c-log 



Qs{Y)J ' 



(6) 



where c is a constant. The AGBS dipole model [l3| inter- 
polates analytically between these two regions through 



N{Y, k) ■ 
with 

Tdii = exp 
and 

ircd = log ( 1 
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Once it describes the evolution for a quark-antiquark 
amplitude, the BK solutions are written in the funda- 
mental representation for the Wilson lines U of the gauge 
fields in the target 
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NF{Y,T)^l-—Tr{u\Q)U{r))^ 



(10) 



The adjoint amplitude for gluons required in the kt- 
factorization formalism can be written, in the large Nc 
limit, as 



NA{Y,r) ^2NF{r,Y) - Nl{r,Y). 



(11) 



In the case of the AGBS dipole model, however, it was 
shown that one can get the adjoint amplitude Na by a 
rescaling of the saturation scale entering the fundamental 
one [l^. In others words, to get Np one can make the 
replacement Ql {CaICf)QI; with Ca/Cf = 9/4, in 
the expression for A^^, in order to make both the hadron 
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production and DIS amplitude compatibles. As showed 
in Ref. [l6j , this is a good approximation of the usual ex- 
pression (fTTj) within the framework of the AGBS model, 
giving a gluon saturation momentum higher than in the 
gluon case. This is what one would expect in the large 
Nc limit where a gluon is seen as two quarks. The result- 
ing forms for both adjoint and fundamental amplitudes 



in momentum space are depicted in Figure 1(a) 



III. 



HADRON PRODUCTION FROM 
itt-FACTORIZATION 



In the regime of very high energies, for y/s ^ Qs > pt, 
the dominant process contributing in the high energy 
cross section is the gluon production via gluon-gluon fu- 
sion and the subsequent fragmentation of the produced 
gluon. When the scattering process occurs between a di- 
lute projectile and a dense target, with large occupation 
number, the cross section for production of a gluon jet 
with transverse momentura qt can be described by the 
fct-factorization formalism [21[ 



dyiPqt 



Cpq^ J 4 

qt + h 



\qt 



(12) 



where xi^2 = {it / \/s)e^'^ are momentum fractions of the 
incoming gluons and Cp = {N^ — i)/2Nc is the Casimir 
for the fundamental representation. 

Even though the model we used is based on the LO 
BK solutions, it incorporates phenomenologically some 
properties of the NLO dipole dynamics, such as the sat- 
uration exponent, that is expected to be reduced from 
A ~ 0.9 in LO models to A ~ 0.2 — 0.3 from the data. 
Thus, concerning the argument of the strong coupling 
constant in , instead of using a fixed value we allow 
for the running of as{Q) within the LO prescription for 
three light quark flavors 



a.(g') = 



127r 



27 log 



(13) 
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^^QCD 



with the maximum value of momentum of the UGDs as 
scale, i.e., we use Q = max{|(7( + kt\/2, \qt — kt\/2}, and 
setting A|cD = 0-05 GeV^. 

In the large Nc limit, the unintegrated dipole gluon 
distribution in either of the two colliding hadrons can be 
related to the dipole scattering amplitude through [12] 
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FIG. L (a) AGBS amplitude in momentum space and (b) the 
respective unintegrated gluon distribution of a proton. 



Fgix^kt) stems for the adjoint dipole scattering ampli- 
tude. As one can see in Eq. ([1]), the BK amplitude in 
momentum space is subject of a slight different Fourier 
transform, but this can be overcome using (22j 

F,(x, k) = l-iylNa{x, k) + S^^\k)), (15) 

where the delta function can be neglected as the uninte- 
grated dipole gluon distribution ([T4|) vanishes for kt — 0. 
Such prescription was successfully used with the AGBS 
model to describe simultaneously the HERA data on the 
proton structure functions and the inclusive hadron yield 
of charged hadrons at RHIC |16l| . The functional form of 



the UGD is shown in Figure 1 1(b) 



IV. RESULTS 



ip{x,kt) 



2Tr^as{kt 



(14) 



where is the number of colors and St is the transverse 
area of the interacting hadrons (nuclei). The quantity 



The transverse momentum distribution of the pro- 
duced charged hadrons is given by. 



di]d^pt ansd J zP- dyd'^qt 



Dh{z ^Pt/qt,y) , (16) 
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where J[r]] is the Jacobian for the transformation of the 
rapidity y and the measured pseudorapidity j] variables 



yiv,Pt) 



■In 



'\J m"^ + p1 cosh^ rj + pt sinh rj 



p1 cosh rj ~ Pt sinh ?y 



(17) 



with TO denoting the hadron mass. In (|16p . Dh{z = 
Pt/kt, fi) stands for the fragmentation function of the pro- 
duced gluon into hadrons, for which we used the LO KKP 
model [231 with the scale ^ = pt oi the hadron and with 
Zmin = 0.05, as demanded from the momentum sum rule. 
To avoid the divergences appearing for small values of pt 
in the /ct-factorization formula ([12]), we made the change 
Pt — )■ \/ m? + Pt, where the hadron mass to has the same 
value as that used in the Jacobian relating the rapidity 
and pseudorapidity variables. 

The non single diffractive cross section cjnsd is, in prin- 
ciple, model dependent and should be taken from soft 
models such as in Refs. [1^ [2^. However, being the 
physical meaning of Unsd the interaction area, we can 
model it, as done for instance in Refs. [13, |3l|. Fol- 
lowing |3l| . we introduce an energy-dependent interac- 



tion radius fov, 



felog(s), such that T^h^^^ mim- 



ics the non single diffractive cross section. We use here 
two alternatives to fix the parameters a and b: to repro- 
duce the results of the model [2^ , that includes both soft 
and hard contributions, which is shown in Figure [2| for 
two different parameters set obtained in such model; and 
to describe the central production of charged hadrons 
dN / drj\ri=Q at different values of the center of mass en- 
ergy. The best agreement of the AGBS model with the 
data implies a = 0.1 and b = 0.198 and the results of this 
choice are shown in Figure [31 We will show our results for 
the charged hadron yield using the first choice for (Tnsd, 
noting that following the KMR model [l^ the central 
rapidity data cannot be described for a broad range of 
y/s in the framework of the AGBS model. Thus, the last 
choice seems to be the best one to fix the normalization 
of our calculation. Being an overall normalization, this 
choice for the non single diffractive cross section allows 
for a K factor of the order of unity, even though one does 
not expect such values for a LO calculation. 

The AGBS model was fitted to HERA data consid- 
ering the small-x region data, x < 0.01, so that quark 
corrections to the CGC formalism could be ignored. In 
proton-nucleus or proton-proton collisions, however, one 
has to consider corrections due to large- a; effects, which 
were introduced in our calculation through the relation 

^{kt,x)^(^—^y(^Y\{kt,xo), (18) 
\1 — xo J \ X y 

where /3 = 4 and Aq can vary from to 0.2 in the proton- 
proton case. 

In Figures [4(a)] and [4(b)] are shown the transverse 
momentum and pseudorapidity distributions of charged 
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FIG. 2. Comparison between the results from the different pa- 
rameters set of the KMR model 28] for the non single diffrac- 
tive cross section a„sd and the prescription a„sd = Trfc; 
where femax = feo -I- C log(s). 
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FIG. 3. Gentral rapidity distribution of hadrons in p{p + p) 
collisions at different energies. The prescription bmax = feo + 
Clog(s) fairly describes the data for the values bo « 0.15 and 
C « 0.2. For comparison, it is also plotted the power law 
curve that fits the data, which was taken from [JuJijl. 



hadrons at LHC for 0.9, 2.36 and 7 TeV p+p collisions, 
which present a good description of the data. Concern- 
ing the K factors we got a better description of the LHC 
data in comparison with the previous work using the hy- 
brid formalism [l6j . This result is expected once the kt- 
formalism is better suited to deal with central rapidity 
data as that measured by the CMS collaboration, while 
the hybrid formalism describes the forward hadron frag- 
mentation region. It is worth to note that such formalism 
needs corrections due to inelastic scatterings of the pro- 
jectile partons on the target field [s^ that should play an 
important role at midrapidity region and at high hadron 
Pt, while in the forward regions, the inelastic piece has 
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a very small contribution to the production cross sec- 
tion [34| . Such corrections, however, are also included 
in the full one loop calculation to the inclusive hadron 
production recently done in [ssl . l36j , and we leave for a 
future study the phenomenological applications of these 
observables within the traveling wave method of QCD. 

In Figure [5] we show the pseudorapidity distribution 
against the recent preliminary data on proton-lead col- 
lisions measured by ALICE Collaboration The 
nuclear saturation scale is modeled through ^ = 

All^Ql p, where A^g = 20 for lead targets. We have also 
checked the prescription of Ref. [s^l , for which Q^,A — 
{AR\/Rl) (with Ra = 1.12^1/3 - 0.86^-1/3 /w), and 
the description is similar. The last one was already 
used with the AGBS model to describe the nuclear ratio 
R{A/B) = BF^/AFi with a good description of 
the shadowing region. 

Looking at Figure [H] one can see that the proton for- 
ward region, for negative pseudorapidity, is well de- 
scribed by the model. This is expected, once the small-x 
(a;2 <?; 1) effects of the nucleus are encoded in this region. 
Regarding the positive rapidity range of the nucleus, it 
should be important to include a better prescription of 
the nuclear geometry. An implementation of the impact 
parameter in the model might improve the description of 
the nuclear data, once it should embedded into the sat- 
uration scale both the large- a; nuclear effects, like EMC 
and Fermi motion, as well as the small- a; nuclear shad- 
owing effects. 

The result presented here for the proton-nucleus col- 
lisions is similar to that of the IP-Sat saturation model 
[39| . mainly in the case of LHC collisions, even though 
we do not use any term to include the impact parameter 
dependence on the saturation scale. Comparing our pre- 
dictions with the data released recently by the ALICE 
collaboration [l| one can sec that the ratio 



dN/df]\^=^2 

dN/dr]\r,=2 



(19) 



is almost the same as the IP-Sat model one, for which 
the ratio is 1.32 against 1.33 for the AGBS model (the 
KLN Hi] and the rcBK ^ models show the ratios 1.38 
and 1.42, respectively.). Thus one can see that concern- 
ing the rapidity evolution, even a LO model as the AGBS 
one can describe quite well the present data, despite the 
large- a; regions we already discussed. The same is not 
true regarding the transverse momentum distribution, 
where the LO dipole models would fatally fail to describe 
the large pt region, once in this region is mandatory the 
inclusion of higher order terms on the dipole evolution, 
or even include a virtuality evolution. One should still 
stress that all these saturation models are not close to the 
value 1.19 of the experiment [l|, maybe signalizing that 
for the LHC data on single inclusive particle spectra a 
evolution is important as well. 
We can also use the central rapidity data of produced 
hadrons to constrain our parametrization for the non sin- 
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FIG. 4. Pt (a) and rj (b) distributions of the LHC charged 
hadron yield for p + p collisions at 0.9, 2.36 and 7 TeV [241 . 
[25( 1. AGBS parameters of the simultaneous fit 16] of AGBS 
to RHIC [li, O and HERA 17j data. The normalization 
used follows Trfe^axi with 6max taken from Figure [2] and the 
K factors are the same for both observables. We note that 
using the normalization from the Figure [3] we could describe 
the data whit a K factor close to one. 



gle diffractive cross section, as show in Figure [3l and thus 
decrease the values of the K factors. Indeed, we checked 
that is possible to describe the data with K factor of 
order one for all the LHC energy spectrum. The RHIC 
data imply a little larger value as consequence of the nor- 
malization showed in Figure [3l that in the range of RHIC 
energy underestimates the data. 



V. DISCUSSION 

We have calculated the inclusive hadron production 
at LHC within the /ct-factorization formalism using a 
saturation based unintegrated gluon distribution. We 
have previously studied in Ref. II6II the same observ- 
able within the hybrid formalism |47{ for the production 
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FIG. 5. Pseudorapidity distribution of the charged hadrons 
produced in p + Pb coUisions at LHC 



cross section, in a global fit that constrained the param- 
eters model to both the HERA data on the proton struc- 
ture function and the inclusive hadron yield at RHIC. 
Comparing both approaches we could observe that, even 
though these formalisms are designed to describe distinct 
kinematical regions of the data - the hybrid one for the 
forward region and the kt factorized at the central ra- 
pidity region, the results are very similar except for the 
values of the K factors normalizing all the uncertainties 
of the model in the multiplicity cross sections. Indeed, 
the K factors gotten in this work are smaller compared 
with the hybrid formalism, as expected once we are deal- 
ing with the central data the CMS experiment delivered 
[24 [25| . Thus a detailed study of the region of validity 
for each formalism is crucial on the study of saturation 
physics at LHC, once such phenomena can be accessed 
in distinct kinematical regions of the detector, depending 
on the observable and kinematical cuts used. 



The large nucleus effects were studied using the re- 
cently measured proton-lead collisions at LHC. In this 
case the AGBS model can describe quite well proton 
fragmentation region where the nuclear small-x effects 
are present. The nucleus fragmentation region lacks the 
implementation of large-a; effects on the nuclear target 
wavefunction. All in all, one can sec that a LO dipole 
evolution can describe the present data on the pseudora- 
pidity distribution. 

Concerning the AGBS model and its limitations we 
note that it describe quite well the data on p + p colli- 
sions, but do not incorporate the large x physics present 
on nucleus fragmentation region, when applied to p + A 
cases. A way to introduce such effects on the model might 
be the study of its whole transverse dependence, what is 
also valid for all the unintegrated gluons (dipolcs) distri- 
butions coming from the small- a; evolution. Being trans- 
verse momentum dependent (TMD) gluon distributions 
functions, it should be important to employ these models 
on semi inclusive observables probing deeply the trans- 
verse momentum on the initial state. Another source of 
information would come from multiparton correlations, 
once on those processes a more complex color structure 
may be probed. Moreover, it is known that the solutions 
of the BK are not suited to describe both the momentum 
fraction x and impact parameter dependencies, once the 
last is intrinsically a non perturbative effect. Thus it 
might be interesting to model such dependencies with 
parameters constrained by diffractive DIS data, so that 
the model could be used to give a better description of 
the transverse plane on nuclear reactions. 
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